% Load data
clc;
close all;
load('activation_energies.mat');
load('max_growthrate_arrays.mat');
load('monod_KM_aggregate_analysis.mat');
load('concentrations_KM.mat');


%% Plot activation vs. log (KM) at 37 °C

figure;
% Glucose
n = length(Ea_glucose);
plot(log(concentrations_KM_glucose(1:n)),Ea_glucose, 'bo-', 'linewidth', 1);
hold on;

% acetate
n = length(Ea_acetate);
plot(log(concentrations_KM_acetate(1:n)),Ea_acetate, 'ko-', 'linewidth', 1);

% succinate
n = length(Ea_succinate);
plot(log(concentrations_KM_succinate(1:n)),Ea_succinate, 'mo-', 'linewidth', 1);

% LB
n = length(Ea_LB);
plot(log(concentrations_KM_LB(1:n)),Ea_LB, 'ro-', 'linewidth', 1);

% NZCasePlus
n = length(Ea_NZCasePlus);
plot(log(concentrations_KM_NZCasePlus(1:n)),Ea_NZCasePlus, 'go-', 'linewidth', 1);

set(gcf, 'Position', [0 0 400 300]);
set(gca, 'fontsize', 20);
box off;
xlabel('Log(Concentration)');
ylabel('Activation energy (k_BT)');